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Abstract 

The production of antibaryons is calculated in a microscopic transport 
approach employing multiple meson fusion reactions according to detailed 
balance relations with respect to baryon-antibaryon annihilation. It is found 
that the abundancies of antiprotons as observed from peripheral to central 
collisions of Pb + Pb at the SPS and Au + Au at the AGS can approximately 
be described on the basis of multiple interactions of 'formed' hadronic states 
which drive the system to chemical equilibrium by flavor exchange or quark 
rearrangement reactions. 
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1 Introduction 



Ever since the first observation of antiproton production in proton-nucleus |l|, 0, |^ 
and nucleus-nucleus collisions [||, ||, H, |] the production mechanism has been 
quite a matter of debate. Especially in nucleus-nucleus collisions at subthreshold 
energies traditional cascade calculations, that employ free A^A^ production and pN 
annihilation cross sections, essentially fail in describing the high cross sections seen 
from 1.5 - 2.1 A-GeV [^, |1^, |Tl|. Thus multiparticle nucleon interactions |T^ 
have been suggested as a possible solution to this problem. On the other hand it has 
been pointed out that the quasi-particle properties of the nucleons and antinucleons 
might be important for the p production process which become more significant with 



increasing nuclear density. Schaffner et al. found in a static thermal relativis- 
tic model based on scalar and vector self energies- assuming kinetic and chemical 
equilibrium - that the p-abundancy might be dramatically enhanced when assuming 
the antiproton self energy in the medium to be given by charge conjugation of the 
nucleon self energy. This assumption implies strong attractive vector self energies 
for the antiprotons which lead to a reduction of the necessary energy to produce 
pp pairs in the medium by binary quasi-particle interactions. On the other hand, 
such self energies will lead to different spectral slopes of protons and antiprotons as 
pointed out in Ref. |]15| . 

First nonequilibrium calculations within a fully relativistic transport model for 
antiproton production - including p annihilation as well as the change of the quasi- 
particle properties in the medium - have been performed in Ref. [n)|. There it was 
found that according to the reduced antinucleon energy in the medium the threshold 
for p-production is shifted to lower energy and the antiproton cross section prior to 
annihilation becomes enhanced e.g. for Si + Si at 2.1 GeV-GeV by approximately 
a factor 70 as compared to a relativistic cascade calculation where no in-medium 
effects are incorporated. Later on, a couple of relativistic transport calculations have 



been performed [|1§, |T^, |T8], [T^ essentially pointing out that all the low energy data 



from proton-nucleus and nucleus-nucleus collisions are compatible with attractive p 
self energies (at normal nuclear matter density po) in the order of -100 to -150 MeV 
p!9| , pO| , pT] . However, it had been stressed at that time that the high antiproton 
yield might also be attributed to mesonic production channels [|^, 0] since pp 
annihilation leads to multi-pion final states with an average abundancy of 5 pions 



25[| , which e.g. might stem from an intermediate state of 2 p-mesons and a pion. 
With new data coming up on antibaryon production from nucleus-nucleus colli- 



sions at the AGS [|6|, [23, |28| and SPS |2|, |0|, 0, |3|, || the A enhancement fac- 
tors seen experimentally were no longer that dramatic as at SIS energies, however, 
traditional cascade calculations employing free production and annihilation cross 
sections again were not able to reproduce the measured abundancies and spectra 
3^, ^ ^ especially for S, S and f2, 0. Here additional collective mech- 



anisms in the entrance channel have been suggested such as color rope formation 
or hot plasma droplet formation [^. In another language this has been ad- 



dressed also as string fusion |^, a precursor phenomenon for the formation of 
a quark-gluon plasma (QGP). 
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The intimate connection of antibaryon abundancies with the possible observa- 
tion of a new state of the strongly interacting hadronic matter, i.e. the quark-gluon 
plasma, has been often discussed since the early suggestion in Ref. that espe- 



cially the enhanced yield of strange antibaryons - approximately in chemical equi- 
librium with the other hadronic states - should be a reliable indicator for a new 
state of matter. In fact, the data on strange baryon and antibaryon production 
from the NA49 and WA97 Collaborations show an approximate chemical equilib- 
rium ^ with an enhancement of the , yield in central Pb + Pb 
collisions (per participant) relative to pBe collisions at the same invariant energy 
per nucleon by a factor ~ 15. As pointed out in Ref. |48] the data on multi-strange 
antibaryons at the SPS seem compatible with a canonical ensemble in chemical 
equilibrium. At AGS energies of 11.6 A-GeV/c, furthermore, a high ratio of A/p of 
3.6til has been reported |^ for central collisions of Au + Au, that is not described 
by any approach so far. Strange flavor exchange reactions help in creating multi- 
strange antibaryons, however, the latter strangeness enhancement factors could not 
be described within traditional transport or cascade simulations without additional 
assumptions such as enhanced string tensions or reduced quark masses P7| . 

In a more recent paper Rapp and Shuryak have taken up again the idea of 
multi-meson production channels for baryon-antibaryon pairs to describe the 
antiproton abundancies in central Pb + Pb collisions at the SPS by introducing 
additionally a finite pion chemical potential which helps in enhancing the multi-pion 
collision rate. Later on, Greiner and Leupold [51| have applied the same concept 
for the A production by a couple of mesons including a or (for the s quark). 
However, such estimates remain schematic unless fully microscopic multi-particle 
calculations support or disprove such suggestions. The problem here is that most of 
the transport models include only binary reactions in the entrance channel whereas 
the final channel of an energetic collision may well consist of many hadrons emerging 
from the decay of strings that are excited in the initial reaction. Thus, as has been 
pointed out quite often |^, 0, detailed balance is not included on the many- 
particle level leading to an improper equilibrium state for large times {t — > oo). 

In this work we will address two separate questions: the first one is of more formal 
nature and related to a transport approach that properly takes into account reactions 
of 2 hadrons — >■ n hadrons and vice versa employing detailed balance (Section 2). 
The second one addresses a suitable covariant scheme for the calculation of such 
multi-particle reactions in transport models. The method and its implementation 



in the hadron-string-dynamics (HSD) transport approach [^, |5^ will be described 
in Section 3. A first application of this novel approach is devoted to the problem 
of antibaryon production in relativistic nucleus-nucleus collisions by multi-meson 
reaction channels. Respective calculations and studies at SPS and AGS energies - 
with a focus on antiproton abundancies - will be presented in Section 4 whereas 
Section 5 concludes this work with a summary. 
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2 Generalized transport equations 



In this Section a brief description of the relativistic transport model is given with 
emphasis on a new development, i.e. the multi-particle reaction dynamics. First we 
summarize (or review) the relevant equations determining the dynamics of baryons 
and mesons and then discuss a flavor rearrangement model for baryon-antibaryon 
annihilation and production in the extended HSD transport approach. 



2.1 Hadron transport and multi-particle transitions 

Since the covariant transport approach for binary reactions has been extensively 



discussed in Refs. and in the reviews |2y, |56|, ^ we only recall the basic 

equations that are relevant for a proper understanding of the results to be reported 
in this study. 

For the discussion of the general collision terms the hadron self energies will be 
discarded for transparency^, such that the transport equation in the cascade limit 
reduces to 

p^d^F,{x,p) = {podt + p- dr)F,ix,p) = Ila, (1) 

where Fi{x,p) is the Lorentz covariant 8 dimensional phase-space distribution func- 
tion for an off-shell hadron with quantum numbers i, i.e. 

F{x,p) = A{x,p)f{x,p), (2) 

where A{x,p) denotes the hadron spectral function while f{x,p) describes the oc- 
cupation probability in phase-space. The off-shell propagation of hadrons leads to 
additional terms in the l.h.s. of (1) that describe the change of the spectral func- 
tion Ai{x,p) during the propagation. These terms are omitted here since the actual 
calculations will be performed in the on-shell limit. For details on the off-shell 
propagation of hadrons in the medium the reader is referred to Refs. [^, |6^ . 



We now turn to a discussion of the collision term I^^n, which includes the new 
elements to be presented below. In the most general case it is a sum of collision 
integrals involving n ^ m reactions. 



I- coll 



^<^oii[n ^ m]. (3) 



The general form for off-shell fermions with spectral functions Ai{x,pi) in case of 
2-body interactions is given by [p9| 



lU^ 2] = 

1 1 /■ 

2 f2 J ^''^'2C?Vc?V Ai{x,p)Aj{x,p2)Akix,p3)Ai{x,p4) 

j k,l V / 

xW2,2iP.P2;t,j\P3,P4;k,l) (27r)^ 6\p''+pi^-p^,-jf,) 



^The actual transport calculations, however, include hadron self energies which optionally can 
be switched-off. 
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^[fk{x,P3)fl{x,p4)il - fi{x,p)){l - fj{x,p2)) 

-fi{x,p)fj{x,P2){l - fk{x,P3)){l - fi{x,P4))]. (4) 

This collision integral describes the change in the 8 dimensional phase-space distri- 
bution function Fi{x,p) due to the collisions of two baryons with momenta and 
P2 and discrete quantum numbers i and j, respectively, whereas the two fermions 
in the final state of the reaction are labeled by their momenta p^ and p^ and dis- 
crete quantum numbers k and I. The 5— function guarantees energy and momentum 
conservation in the individual collision while W2,2{p,P2',hj \ P35P4; denotes the 
transition probability (or transition amplitude squared) for this reaction which in 
case of fermions includes the antisymmetrization. 

The generalization of (§) to n m interactions is straight forward and reads: 

^^rEE ((^) ^ / {Ud% M^,P^?j [f[^d%A,ix,p,)^ 

n m 

xAix,p) Wn,m{P.Pf.h^\Pk-A) {27lf 5\p^ + P' ' Pk) 

j=2 k=l 

m n n m 

^[fi{x,p) n fk{x,Pk) n fji^^Pj) - fi{x,p) n fA^^Pj) n fk{x,Pk)\- (5) 

fc=l j=2 3=2 k=l 

In Eq. the quantities / denote Pauli-blocking or Bose-enhancement factors as 

f = l + Vf (6) 

with 7] = 1 for bosons and rj = —1 for fermions, respectively. The indices u and 
A stand for the set of discrete quantum numbers in the initial (except for particle 
i) and final states, respectively, and denotes a statistical factor that takes into 
account the number of identical fermions and bosons in the initial and final states. 
In the following we will assume that the transition probabilities Wn^m are evaluated 
with respect to asymptotically free antisymmetrized fermion many-body states and 
symmetrized meson states such that N'^ = 1. 

In the on-shell quasi-particle limit the spectral functions, furthermore, are given 

byQ 

A{x,p)=2n6ip'-M^) (7) 
which leads to the 2-body collision integral for particle i (in case of fermions) as: 



1 1 f d^P2 d?Pz d?Pi 

. (2^ J 2E^2E^2E^ 



IIou[2^2]--YYt^, j 



xW2,2ip,P2;i,j\P3,P4;k,l) (27r)4 S^pf^+pi^-p^^-pl) 

^[fk{x,P3)fl{x,P4){l - fi{x,p)){l - fj{x,p2)) 
-Mx,p)fj{x,P2){l - fk{x,P3)){l - fl{x,p4))], 



^In general the spectral functions for fermions differ from that of bosons |59|; here baryons and 
antibaryons with different spins are treated as independent Klein-Gordon particles. 
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where the energy p\ is rewritten as Ek. Note that the po = Ei integration over 
the spectral function Ai{x,p) appears on both sides of the transport equation and 
cancels out in the limit Fj 0. The on-shell version of the collision integral 
then reads 



n m 

xW^,m{P^Pf 2, u I pu- A) (27r)4 5\p^ + E^i - E Pk) 

3=2 k=l 

m n n m 

[fi{x,p) n fk{x,Pk) n fji^^Pj) - fii^^p) n fji^^Pj) n fk{x,Pk)]- (9) 

k=l j=2 j=2 k=l 

For large times {t —>■ oo) all collision integrals vanish, which implies that 'gain' and 
'loss' terms become equal in magnitude. 

The number of reactions in the covariant 4-volume d^rdt = dVdt is obtained by 
dividing the gain and loss terms in the collision integrals by the energy p^ = Ei, 
integrating over d^p/{2'Kf and summing over the discrete quantum numbers i. For 
the case of fermion two-body collisions this gives (using p = pi) for the 'loss' term 



dNcoii[2 2] _ y^y^ 1 r (^jh (fp2 d^ d^pi 
dtdV (27r)i2y 2Ei 2^2 2E^ 2E^ 

y<W2,2ipi,P2;iJ \p3,P4;k,l) (27r)^ S^iPi+P2-P3-Pi) 

^[fiix,Pi)fjix,P2)il - fkix,P3)){l- fi{x,P4))]. (10) 

In case of n ^ m processes this leads to 

dN^ou[n m] ^ / 1 \ "^"^ f -pV ^ yr <^Pk 

dtdV " ir V (2^) V J M l}^ 2E, 

n m n m 

xTy„,„(p,;z,z/ I pk-A) (2vr)^ 5\Y.P' " EPfe)(n fA^^Pi) 11 h{x,Pkmi) 

j=l k=l j=l k=l 

and in case of m — ^ n processes to 

dNcou[m n] ^ v-v- / 1 Y^"" f yr <^Pi -TT (^P^ 
dtdV " tr Y V (2vr)3 j J fA ^E, i\ 2E, 

n m m n 

xWnAPr.h y I Pk\ A) (27r)^ b'i^^p] - E^Ddl f^ix^p^) n /,(x,p,)).(12) 

3=\ k=l k=l j=l 

For the phase-space configurations of interest in this study the Pauli-blocking or 
Bose-enhancement terms fk are ~ 1, which implies to replace the quantum statistical 
ensembles by classical ones. In this limit the integrals over the final momenta can 
be carried out provided that the transition probabilities Wn^m do not sensitively 
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depend on the final momenta pk- Employing the definition of the n — body phase- 
space integrals for total 4- momentum [|6l|l, 



(13) 



one obtains the recursion relation pT| 



i^„(P^ Ml, .., M„) = -i^ J ^P„_i(P'^ - p^; Ml, .., M„_i 



(27r 



(14) 



Note, that the phase-space integrals are of dimension GeV^" or (1/fm)^" ^. In- 
serting (p!3D this gives in case of n — m processes 



dNcoii[n m] 



dtdV 



EE 



n 



1 • 1 2Pt 



A v(2vr)V 

xPV„,^(P) P„(P^ = Y.pf,M^,..,M^) n /,(x,p,) 



(n m)J^ I n /i(a;,Pi) 



(15) 



where Mi,..,Mm stand for the masses in the final state, and in case of m — >■ n 
processes 



dNcoii[m n] 
dtdV 



EE 

i,^ A 



1 



f27r)3 



n 

\k=l 



d^Pk 
2Ek 



k=l 



Kk=l 



= EE((^)"/ (n/^ P{m-^n)r [flh{x,Pk)y (16) 

For fixed sets of quantum numbers {i, v) and A in the initial and final states this 
relates the integrands P{n —>■ m) in ([l5|) and P(m — > n) in (|T6|) for individual 
scatterings as (dropping the indices for quantum numbers): 



P(m n) 
Pin — > m) 



m 1 

n 

L/c=l 



2Pfc 



i?n(P^ = ELiPfc;Mi,-,M„) 
Rm{P'' = J:T=iPj;Mu..,Mj 



(17) 



if essentially depends only on the invariant energy = vP^. Note, that the 
r.h.s. ' of (|13) is in units of GeV^^"""^) or fm3("'-") such that a factor (c/V")""™ is 
needed to interpret the quantities as relative 'probabilities'. 

Thus, once the transition probabilities Wn,m are known as a function of ^/s for 
a given set of quantum numbers, the integrand P{n ^ m) in (|15D is determined by 
phase space and the backward reactions in ([T6|) are fixed by Eq. ([T7|) . 
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2.2 Antibaryon annihilation and recreation 

In the following the processes BB ^ m mesons are discussed, which are of rel- 
evance for annihilation of antibaryons on baryons and the recreation of BB pairs 
by m meson interactions. The 4-differential collision rate for baryon-antibaryon 
annihilation (1 + 2— > 3,.., m + 2) then is given by 

dNcoii [BB —>■ m mesons] / 1 V / ^^P^ ^^P"^ 

The integrand is related to the annihilation cross section cTann.iV^) ^'^^ ^ baryon- 
antibaryon pair with quantum numbers i,j as 



W2,m{Pl+P2;i,j;>^m) Rm{Pl + P2', . . , Mm+2) 



(v^) = 4E^E2 (v^) (19) 



with the Lorentz-invariant relative velocity pi], 52 



^'X{s,MlMl) 

= 2E,E2 ' 

involving 

Kx,y,z) = {x - y - - Ayz. (21) 

In (|l^) the sum runs over the final meson multiplicity (m ^ 2, ... , 9) in the final 
state and over Am which denotes all discrete quantum numbers of the final mesons 
for given multiplicity m. 

Note, that by summing (p!8| ) additionally over m, but keeping the quantum num- 
bers i,j fixed, one arrives at 



dN. 



coll 



^ J d^Pid^P2 Vrel{pi,P2) (^anniVs) fi{x,Pi)fj{x,P2). 



(22) 



dtdV (27r) 

If the product of the relative velocity and the cross section, i.e. VreiO'ann, is approxi- 
mately constant (see below) the integrals over the momenta in (^2]) give the classical 
Boltzmann limit 

dN'''^ 

—j^ = < Vrel (Tann > pi{x)pj{x), (23) 

dtdV 

where Pi{x) is the density of the hadron with quantum numbers i. 

The number of reactions per volume and time for the back processes is then 
given by (A^ = ki,..,km) 

dNcoii[ni mesons BB] / 1 A™" /" /"Ji^ d^pA 



dtdV |J5i(2^)V ' [Us ^E,J 

m+2 /m+2 \ 

xl^2,™(v^;^,J, Am) R2{P^ = E Pk;Mi,M2) n /^(^'Pfc) > (24) 

fc=3 \fc=3 / 
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assuming W2,m{y/s; i, j, \m) to depend only on the available energy y/s and conserved 
quantum numbers. 

To proceed further, some simplifying assumptions have to be invoked to lead to a 
tractable problem for antibaryon annihilation and production. Experimentally, the 
differential multiplicity in the pions from pp annihilation at low -^i above threshold 
can be described as 

1 CA^ - < >)2 



with an average pion multiplicity of < N^^ >~ 5 and D = 0.95 p5[ . 

This observation is reminiscent of flavor rearrangement processes in the BB 
annihilation reaction to vector mesons and pseudoscalar mesons, e.g. p + p + vr or 
uj + u + TT^, where the p and u 'later' decay to 2 or 3 pions, respectively. In this 
picture the p + p + vr final channel in pp annihilation is the dominant process leading 
finally to 5 pions. Alternatively, the u + uj + channel leads to 7 pions in the final 
channel which will appear on the scale of the a;-meson lifetime. Three pions are 
obtained in the direct 3 pion decay which, however, is substantially suppressed at 
higher ^/s due to spin multiplicities (see below). 

For the problem of interest we thus employ a quark rearrangement model for 
BB annihilation to 3 mesons as illustrated in Fig. |l|, where the final mesons 
may be pseudoscalar or vector mesons, i.e (vr, K, rj) or (p, u, K*, 0), respectively. In 
this model there is no creation of an ss pair for nucleon-antinucleon annihilation 
due to the conservation of constituent quark flavors. As can be extracted from the 
detailed experimental cross sections given in Ref . , such processes are suppressed 
by more than an order of magnitude. In principle, such channels can additionally 
be included in the model, however, the backward reactions then are also suppressed 
by the same relative factor according to detailed balance. For the study of interest 
such 'small' channels are discarded. 

In the following, the quantum numbers denoted by will be separated into 
different channels c, that can be distinguished by their mass decomposition, and 
degenerate quantum numbers such as spin multiplicities and isospin projections. 
In the latter sense the sum over the final quantum numbers in ([T8|), (2^) then 



includes a sum over the mass partitions c = (M3, M4, M5) , a sum over the spins of the 
mesons and a sum over all isospin quantum numbers, that are compatible with charge 
conservation in the transition. The probability for a channel c = (M3, M4, M5) then 
reads 

Pe(v^; M3, M4, Ms) = NsiV^s) RsiV^s; M3, M4, M5) N'j,^, (26) 
where the number of 'equivalent' final states in the channel c is given by 

Nf^n = (2^3 + ms, + 1)(2.5 + 1)§^. (27) 

In ( p7D Sj denote the spins of the final mesons, Fiso is the number of isospin pro- 
jections compatible with charge conservation while Ni^, is the number of identical 
mesons in the final channel (e.g. = 3 for the tttttt final channel). This com- 
binatorial problem for the final number of states A^/m is of finite dimension and 
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easily tractable numerically. For each mass partition c = (M3, M4, M5) the decay 
probability then is given by the 3-body phase space and the allowed number of final 
states N'^in since the absolute normahzation - described by N^{^) - is fixed by the 
constraint I]c-Pc = 1- 

As an example let us consider the problem of nucleon-antinucleon annihilation, 
where the following final meson channels contribute, 

(1) TTTTTT (2) VTTrp (3) TXTIU (4) TTpp (5) H pUJ (6) HUJUJ , (28) 

excluding 3 vector mesons in the final channel. According to (|26|) the distribution 
in the final number of pions (including the explicit vector meson decays to pions) 
can be evaluated as a function of ^/s since it only depends on the phase space and 
the number of possible final states A^jj„ in each channel c. The numerical results are 
displayed in Fig. |^ for 2.3 GeV < ^/s < 4 GeV in comparison to the parametrization 
(p5|) (solid line). Here the horizontal bars indicate the range of A^jr-pion probabilities 
when varying ^/s from 2.3 to 4 GeV. Obviously, the simple phase-space model (p6D is 
in a fair agreement with the experimental observation. For related or more extended 



models for pp annihilation the reader is referred to Ref. |25 . 

For the backward reactions, i.e. the 3 meson fusion to a BB pair, the quarks and 
antiquarks are redistributed in a baryon and antibaryon, respectively, incorporating 
the baryons A^, A, A, S, S*, S, S*, Vt as well as their antiparticles. In line with (p6[) 
the relative population of states (with the same quark content) is determined by 
phase space, i.e. 

P,,{^s-M^,M2) = N2{^) /?2(v^;c'= (Mi,M2)) (2si + l)(2s2 + l) = 

A^2(v^) i?2(v^;Mi,M2) Ni, (29) 

where A^^' now denotes the number of final states for the particular mass channel c' 
in the backward reaction. The absolute normalization N2{^/s) is fixed again by the 
constraint J2c' Pc' = 1- 

As an example consider the reactions tt^tt+tt^ or tt^p+tt^ or n^p^p^ (and isospin 
combinations), i.e. ud + du + ud {uud) + {udd): here the final states may be either 
p + n, A~ + n, p + A° or A~ + A'' within the Fock space considered. Note, that the 
final states with a A- resonance are favored due to the spin factors in (^91), however, 
somewhat suppressed by the 2-body phase-space integral R2{\/s) for low a/s. 

One is thus left with the BB annihilation problem 

dNcou\BB ^ mesons] 1 f d^Pi d^P2 , r\ 



dtdV V ^ (27r)6 J 2Ei 2E2 

xN,{y^) R,{p,+p2;c={M,,M^,M,)) Ar;.„ /,(x,pi)/,(x,p2), (30) 

where (Mi,M2) denote the baryon and antibaryon masses in the channel c' and 
(M3, M4, M5) the final meson masses in the channel c. Eq. (|30|) can be rewritten as 

dNcoii[BB — > 3 mesons] 

dtdV " 

EE / d'pi d'p2 P^W^ /.(x,pi)/,(x,p2) (31) 



10 



with the channel probabihties 

P^fiV^s) = W'^V^s) N.iyTs) R^^p^ + p^- (M3, M4, M5)) iV;,„. (32) 

Note, that by construction we have 

E Pcc'iV~s) = W'^'iV~s) = Vrel Cr.„„(v^)e', (33) 

where v^ei denotes the relative velocity (|20|) and (j ann{.\fs) d is the total annihilation 
cross section for BB pairs of channel d . 

The backward invariant collision rate is given by 

dNcoii [3 mesons BB] -^-^ 1 f f tt '^^Pk ' 



xiV2(v^) i?2(Epfe;c' = (Mi,M2)) iV^' fn fk{x,Pk)] . (34) 

A:=3 \fc=3 / 

Using Eq. (|T^, the relation ([T^ for 3 mesons in the final state and (^) one arrives 

at 

dNcou[3 mesons -» BB] _ „ „ 1 r d^P3 d^p^ d^p^ ^ 
dtdV ~^^{27t)^J 2^3 2^4 2E5 ^ ^ 

( r.\ ^2(yi) R2{P^:c' = {M,,M^)) Nj 
xv^ei aws), ^^^^^ ^ ^ ^^^^ ^^^^ 

^h{.x,p3)fi{x,pi)h{x,p^) (35) 

for the backward reaction 3 + 4 + 5^1 + 2. Eq. ( |35D can now be rewritten as 

dNcoii^ mesons BB] 

dtdV ~ 

EE j ^^P3 d^pi d^p5 P^diVs) f3{x,P3)h{x,P4)h{x,P5) (36) 
with the 'transition integrand' 

^^'^ - 2E^"-' iV3(v/^)i?3(P;c=(M3,M4,M5)) iV^' ^^^^ 

which is of dimension GeV~^ or fm^. 



3 Numerical implementation 

For a reformulation of the 'transition integrands' (specified in (P7[)) in a test-particle 
representation one has to recall that the average density of a meson with quantum 
numbers k is obtained by integration over momentum as: 

1 ^.3. 



M^) = j^;^ I d p fkix,p), (3J 
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where e.g. charge, strange flavor content, total spin and spin projection are specified 
by the discrete quantum number k. The conversion formula thus reads: 



(2vr) 



j d?p Ux^p) ^ ^ , (39) 



i e dV 



where dV is a (small) finite volume and the sum runs over all test particles in 
the volume dV with quantum numbers k. The number of BB annihilations in the 
volume dV during the time dt is thus given by 



NbB = -^ J2 Vrel{hj)(^anniVSij) (40) 
i,j e dV 

with the invariant energy squared 

Si,j = {Pl+P2)\ (41) 

where pi,p2 denote the 4-momenta of the colliding BB pair. The relative velocity 
Vrei{i,j) is given by (pO|) while the annihilation cross section cfanniVs), furthermore, 
has to be specified for all baryon-antibaryon pairs. This cross section is rather well 
known for nucleon-antinucleon reactions p5|, however, the channels involving 

A, E, S*, H, H*, Q~ baryons or their antiparticles are not available experimentally by 
now and have to be modeled to some extent. 

For guidance we recall that the product Vrei crann{^/s) ~ 50 mb for a wide range 
of energies ^/s in case of pp annihilation. This is demonstrated explicitly in Fig. |^, 
where the experimental annihilation cross section for pp from Ref. is compared 
to the approximation 

C^ann(v^) = (42) 
Vrel 

which holds well in the dynamical range of interest. We thus can adopt the Boltz- 
mann limit ( P5| ) to estimate the BB annihilation time at nucleon density p as 

Tann. ~ (S/m^)"' ^ 1-2 ^ [fm/c]. (43) 

p 

In this work we will proceed with dynamical calculations in the strangeness sector 
S=0, thus essentially addressing the p abundancies in relativistic nucleus-nucleus 
collisions. In this case the approach formulated above does not involve any new 



parameter or cross section; it is just an extension of the HSD approach [g^, |5J, |65 
to include the BB 3 meson reactions by detailed balance. 

The number of backward reactions by 3 mesons in the test-particle picture in 
the volume dV and time dt according to (|35|) for a given mass channel c' is given by 



iVs™ = ^^^^E^^^|f|;--(1^2).(v/i). 
^, N2{V^) R2{V~s-c' = {M.^M^)) Nj ^ ^^^^^ ^^^^ 



iV3(v^) R^{^s-c={M,,M,,Mu))N' 



fin ij^k e dV 
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where the channel c is defined by the colliding mesons (cf. ([2^)) and the outgoing 
channel c' by the BB pair with masses Mi and M2 and energies Ei and E2, re- 
spectively. In (Q) the summation over the mesons in the volume dV is restricted 
to 2 < j < A; in case of 3 identical mesons (e.g. 3 pions) and to i < j in case of 2 
identical mesons i,j in order to account for the statistical factor A^j^! in Eq. (|27|) . 

Eqs. ( plD and (^) are well suited for a Monte Carlo decision problem, i.e. a 
transition is accepted if the probability P^k is larger than some random number in 
the interval [0,1]. One has to assure only, that all Pijk are smaller than 1, which 
- for a fixed volume dV - can easily be achieved by adjusting the time-step dt. 
This evaluation of scattering probabilities is Lorentz-invariant and does not suffer 
from geometrical collision criteria as in the standard approaches |6^, that imply 
a different sequence of collisions when changing the reference frame by a Lorentz 



transformation For 2^2 transitions it has first been employed and tested by 
Lang et al. in Ref. ||6^; this method is also implemented in the HSD approach, 
where it can be used optionally instead of the standard geometrical collision criteria 



as described e.g. in Refs. [p7| , |66| , |69| . The present implementation in this respect is 
a straight forward generalization of the concept in Ref. to 2 3 reactions. It is 
worth to point out that this numerical implementation is a promising way to treat 
n ^ m transitions in transport theories without violating covariance or causality. 
In case of infinitesimal volumes dV and time steps dt it gives the correct solution to 
the many-particle Boltzmann equation. 

For the actual numerical calculations a dynamical time-step size dt is employed 
which on average amounts to dt ^ 0.5/7cm [f^i/c], where is the Lorentz factor 
in the nucleus-nucleus cms, i.e. ~ 9.3 for collisions of Pb + Pb at 160 A-GeV. 
The volume dV is chosen to be dV = A dz with the transverse area A = 
and dz = S/'^cm [iTa]. Variations of these parameters within a factor of 3 do not 
change the numerical results to be presented below. In order to avoid numerical 
artefacts, that are due to the finite volume dV, a BB pair, that has been produced 
by meson fusion, is not allowed to annihilate on each other again without performing 
an additional collision in between. On the other hand, the mesons stemming from 
a particular BB annihilation are not allowed to fuse again with the same partners 
for the backward reaction, if no intermediate extra collision has occured. 

As a numerical test the number of collisions in a single box of volume 10 fm^ 
during the time dt =1 fm/c has been calculated with spatially uniform phase-space 
distributions given by a classical system of hadrons in thermal and chemical equi- 
librium, i.e. 

Mp) = ^^^^^^^^M-E,{P)/T) (45) 

with s and / denoting spin and isospin, respectively. The particles taken into account 
are A^, A and their antiparticles and vr, p, uj on the meson side in the strangeness 
sector S=Q. The numerical results for the number of BB annihilation collisions 
(— *■ vrpp) are shown in Fig. ^ in terms of the dashed line as a function of ^/s, which 
corresponds to the invariant energy in an individual collision. As can be seen from 
Fig. ^ the dashed line very well coincides with the solid line that corresponds to the 
energy differential number of Tipp collisions for the backward reactions. Thus the 
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numerical scheme employed well reproduces the detailed balance relation in thermal 
equilibrium for a given channel combination cc'. Without explicit representation we 
mention that the detailed balance relation is fulfilled for all channel combinations 
cc' specified above. 



4 Nucleus-nucleus collisions 
4.1 SPS energies 

The most complete set of data on antibaryon production in nucleus-nucleus collisions 
is available from the NA44 J7|, NA49 and WA97 |3l|, |32| Collaborations for 



Pb + Pb collisions at SPS energies of 160 A-GeV, that allow for stringent tests 
of the dynamics proposed. Since the extension of the HSD transport approach is 
described in the previous Section and no new parameters enter into the calculations, 
we proceed with the actual results. 

4.1.1 Cascade calculations 

Though antiproton self energies have been found to be important in nucleus-nucleus 



collisions at subthreshold energies |19, 20, 21 1, we start with cascade calculations 



because the initial invariant energy per nucleon is large compared to the BB thresh- 
old. However, before coming to the antibaryon abundancies the performance of the 
HSD transport approach has to be tested in comparison to experimental data for 
baryons and mesons. Since detailed differential spectra for protons, hyperons, pi- 
ous and kaons have been presented in Ref. |^ in comparison to the experimental 
data for central collisions of Pb + Pb at 160 A-GeV, we concentrate here on parti- 
cle abundancies as a function of 'centrality'. The latter is defined by the number 
of participants Apart, which is extracted from the transport calculation at impact 
parameter b as 

Apartib)=2A-Noib), (46) 

where NQ{b) is the number of nucleons that were not involved in any hard scattering 
process. 

The average number of charged pions < vr >= (< 7r+ -|- vr" >)/2 (divided by 
Apart) is shown in Fig. 5 as a function of Apart in comparison to the data from 
Ref. for Pb + Pb at 160 A-GeV. The number of pions per participant nucleon is 
found to be approximately constant within 10 % as a function of centrality; there is a 
slight trend in the data as well as in the HSD calculations (solid line) for an increase 
of < vr > /Apart for central collisions. However, the HSD transport calculations 
overestimate the pion abundancy by about 20 %. Such an overprediction of pion 
multiplicities in heavy systems occurs in other transport approaches as well 



- or is even higher - and is not well understood so far. At SIS energies of 1-2 
A-GeV it has been argued in Ref. that a quenching of nucleon resonances at 



baryon densities p > po might be responsible for the relative pion suppression seen 
experimentally in heavy systems, but it is not yet clear if such a mechanism will 
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also explain the discrepancies at SPS energies. We thus have to keep in mind this 
overprediction of pions especially when comparing particle ratios (see below). 

The effect of antibaryon annihilation and reformation by means of meson fusion 
channels is shown in Fig. 6, where the < p > / < it > ratio is displayed for Pb + Pb 
at 160 A-GeV as a function of Apart- The dashed line shows a calculation without 
including annihilation channels in the transport calculation, the dotted line gives the 
results including the annihilation to mesons while the solid line is obtained when 
including both, annihilation and meson fusion channels by detailed balance. The 
first point in Fig. 6 gives the numerical result for A^A^ interactions at Tiab = 160 GeV 
{Apart = 2). As shown in Table 1 of Ref. [^, the description of the HSD approach 
for pp interactions is well in line with the data on 7r+, vr", vr", AT"*", K~, K^, A + S'', 
A + S° as well as p and p multiplicities at SPS energies, such that this point may 
serve for reference to the particle abundancies in the elementary A^A^ interaction. 
As seen from Fig. 6 the < p > / < tt > ratio slightly drops with centrality for 
peripheral reactions, however, stays approximately constant for Apart > 150 when 
neglecting annihilation. Thus antiprotons in central Pb + Pb collisions are produced 
less frequent than pions by about 30 % relative to the A^A^ interaction in vacuum. 
The effect of annihilation (dotted line) sets in already for peripheral reactions and 
amounts to a factor ~6-7 suppression for central collisions. The latter suppression 
is sensitive to the formation time tf of the antibaryons and becomes larger (smaller) 
for shorter (longer) Tp. In the actual calculations we have used Tp = 0.8 fm/c for 
all hadrons When including annihilation as well as meson fusion channels by 
detailed balance, the < p > / < n > ratio (solid line) becomes again close to the 
calculation that does not include annihilation nor meson fusion reactions to BB 
(dashed line). Thus on average the annihilation channels are almost compensated 
by the reproduction channels. 

In order to get some idea about the dynamical origin of the approximately con- 
stant p/tt ratio shown in Fig. ^ by the solid line the reaction rate B + B ^ mesons 
(dashed histogram in Fig. ^ ) is compared to the backward reaction rate (solid 
histogram in Fig. |^ for a central collision of Pb + Pb at 160 A-GeV. Fig. |^ demon- 
strates that both rates are comparable within the statistics. Thus an approximate 
local chemical equilibrium is established very fast between the nonstrange antibaryon 
degrees of freedom and the nonstrange mesons vr, p and u. The latter fact is sup- 
ported by the absolute number of reactions BB mesons which is about 4-5 times 
higher than the final number of antibaryons per event. On the other hand, the num- 
ber of backward reactions is ~ 96% of the number of annihilation reactions leading 
to a small net absorption of the antibaryons produced initially by baryon-baryon 
(~ 73%) or meson-baryon (~ 27%) inelastic collisions. Only ~ 12% of the final an- 
tibaryons stem from 'hard' baryon-baryon or meson-baryon reactions; the dominant 
amount of final p's (~ 88%) are from the 3 meson fusion reactions indicating that 
the memory with respect to the initial 'hard' collision phase is practically lost. 

One might worry about the sensitivity of these results to the annihilation cross 
section crann{y/s) that so far has been taken to be the free cross section (in the 
parametrization from Ref. |T^). However, this quantity might change in the medium 
due to screening effects. It is clear that any enhancement of this cross section will 
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lead to an even faster equilibration in the light flavor degrees of freedom and to a 
more perfect chemical equlibrium. Thus numerical calculations have been performed 
for central Pb + Pb collisions at 160 A-GeV by assuming that aann is reduced by a 
factor of 2. This leads to a reduction of the total number of annihilation reactions 
(and backward reactions) by a factor of 2 - which essentially reduces the numerical 
statistics - but leaves the conclusions unchanged. Thus the flndings from Figs, ^and 
1^ are robust against 'reasonable' modiflcations of the in-medium transition rates. 

The numerical abundancies for mesons and antiprotons relative to the aver- 
age charged pion multiplicity < vr >=< 7r+ + vr^ > /2 are displayed in Fig. || for 
Pb + Pb at 160 A-GeV as a function of Apart- The comparison of the calculations 
with the data from Ref. [0 indicates that the dependence of the multiplicities 
on Apart is roughly met - except a single point for K~ at rather peripheral reactions 



- in line with the earlier analysis in Refs. which concentrated on central 

collisions in this system. However, the p/vr ratio is lower by about a factor of 2 com- 
pared to the data of the NA49 Collaboration that also include the feeddown from 
A and S° due to the weak interaction. Since the latter contribution is presently 
unknown one might either speculate that the A + I]° abundancy is comparable to 
the antiproton abundancy or that antiproton self energy effects might be responsible 
for the experimental observation. 

4.1.2 Antiproton self energies 



As shown in Refs. |T6| , ^ the production of antiprotons at SIS energies of 

1.4-2.1 A-GeV as well as in p + A reactions is described by adopting attractive self 
energies for the antiprotons in the range of -100 to -150 MeV at normal nuclear mat- 
ter density po- Especially in p + A reactions the backward production channels by a 
couple of mesons are statistically irrelevant as can be easily checked by the transport 
model described above. The question thus arises if such 'established' antiproton self 
energies for densities 1-3 po might also be responsible for an enhancement of the p 
yield at SPS energies in the Pb + Pb system. 

To examine this possibility we show in Fig. ^ the time evolution of the baryon- 
density in a central cylinder of transverse radius Rt = 5 fm, that has moving bound- 
aries with the expanding hadronic system in longitudinal direction. The solid line 
denotes the average density of 'formed' baryonic states whereas the dashed line rep- 
resents the net quark density pq/3 which merges with the baryon density in the later 
expansion phase. Both densities have been evaluated in the cms rapidity interval 
l^ylcm < 1 in order to exclude spectator nucleons and to gate on midrapidity physics. 
The difference between the solid and dashed line in Fig. ^ has to be attributed to 
'non-hadronic' states, which in the HSD approach are quarks and diquarks (as well 
as their antiparticles) that constitute the ends of 'strings' or continuum excitations 
of the hadrons. As can be seen from Fig. ^ the 'non-hadronic' phase lasts about 2.2 
fm/c which is roughly the diameter of the target (2i?) divided by the Lorentz-factor 
^ 9.3 plus the hadron formation time Tp = 0.8 fm/c 



2R , , 

Tnonhad. ~ V Tp- (47) 



16 



Then a mixed phase of 'partons' and 'formed' hadrons comes up which practically 
ends around 6 fm/c where all continuum excitations have merged to hadrons or 
hadronic resonant states. Quite remarkably, the density of 'formed' baryons is about 
2po at the beginning of the pure hadronic expansion phase. Now the BB annihilation 
rate (in Fig. |^) starts with a maximum around 3 fm/c - when the 'baryons' also 
start to hadronize - along with the meson fusion reactions that last up to about 8 
fm/c. The characteristic time scale for the decrease of the annihilation/fusion rate 
is Tprod ~1.6 fm/c for central Pb + Pb reactions at the SPS as extracted from Fig. 
1^ using an exponential ansatz. During the BB production phase by meson fusion 
from 3-8 fm/c the density of baryons drops from ~ 2.5po to po- Employing (43) 



the time scale for annihilation at these densities changes from 0.5 - 1.2 fm/c, which 
is considerably shorter than the characteristic production time scale of 1.6 fm/c. 
The approximate chemical equilibration between mesons and antibaryons thus is 
no surprise according to these simple 'classical' estimates. On the other hand, the 
densities of 2.5 po ^ Po are very similar to the densities probed in heavy- ion reactions 



at SIS energies from 1-2 A-GeV |21, 57 



In order to investigate if p self energies might be responsible for the experimental 
antiproton yield, a scalar attractive p potential of the form 

Up{p) = -a-^ (48) 
Po 

is assumed with a = 100 MeV, where p denotes the density of 'formed' baryons. 
This amounts to produce and propagate antiprotons with a density-dependent mass 
M* = Mq + U{p). For more technical details the reader is referred to Refs. [1^, [19[. 
Note, that when introducing self energies for particles by averaging the latter over 
many events one no longer performs microcanonical simulations since the field energy 
is shared between different events. However, quark flavor conservation still holds 
exactly in each event such that the physical system corresponds to a canonical 
ensemble. 

The numerical results for the p/ < it > ratio in Pb + Pb collisions at 160 
A-GeV are displayed in Fig. [l^ by the solid line in comparison to the cascade 
calculation (dashed line) and the experimental data from NA49 (full triangles). 
It is seen that with increasing centrality or Apart there is a slight enhancement of 
the p yield for the potential (^), which vanishes for very peripheral reactions where 
the average baryon density is very small. The experimental data, however, are still 
underestimated significantly. This also holds true when accounting for the 20% 
overprediction of pions in the transport approach (cf. Fig. 5). 



4.1.3 Extrapolation for antihyperon production 

In view of the approximate chemical equilibrium achieved for the u, d quark sector 
between mesons and antibaryons - even for reduced annihilation cross sections - 
one may proceed with speculations on the strangeness (ss-quark) sector, where no 
experimental data on the annihilation cross sections are available. However, in case 
of similar transition matrix elements squared (~ 5 fm^) the YN and YN channels. 
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where Y stands for the hyperons (A, S, S*), will achieve chemical equilibrium with 
the 3 meson system as before, however, with a vr or p, exchanged by a or 
K*,(j), respectively (cf. Ref. |5ll] for the case of 5 meson reaction channels). A 



further step then consists in replacing another pion or p, by a ii" or K* , which 
will bring the S, S* and S, S* in chemical equilibrium with the SU{i) flavor meson 
system. Moreover, in case of flavor rearrangement reactions with 3 strange mesons 
K*, (j) (or any combinations of those) the Q, Cl system might also achieve chemical 
equilibrium. This will imply for the antibaryon to baryon ratios: 



p^K+ \K+ J S ^ \K+ J Q' 



(49) 



The relations ( ^91) are easily obtained for a grand canonical ensemble at 'high' tem- 
perature T where Fermi and Bose distributions coincide with the Boltzmann distri- 
bution. In this limit the ratios of particles to antiparticles (apart from the tempera- 
ture T) only depend on the chemical potential p^ for light quarks and /ig for strange 
quarks, i.e. 

P _ exp(-(£'p + 3/i,)/T) 



p exp-{Ep-3fig)/T) 

K_ ^ exp(-(£;j^ -p, + pg)/T) 
K exp{-{EK + ^s- Pq)/T) 

A ^ exp(-(EA + /i. + 2/ig)/T) _ 
A exp(-(EA - p, - 2/ig)/r) 

S ^ exp(-(Eg + 2/i, + /z,)/r) _ 
S exp(-(EH - - p5)/T) 

n exp(-(Ef^ + 3/i,)/T) 



= exp(-6/ig/r), 
= exp(2(^, -pq)/r), 
exp(-(2p, + 4;U,)/T) 
exp(-(4p, + 2/i,)/T) 
= exp(-6p,/T), (50) 



n exp{-{En-3fis)/T) 

where Ex, Ex denote the energies of particles and antiparticles, respectively, that are 
the same when neglecting self energies. Note, that the relations (^) also result from 
the quark condensation model in Refs. [|73| , [f^ , however, involve a slightly different 



physical picture. In the latter approach the mesons, baryons and antibaryons emerge 
from an equilibrium QGP state by condensation under the constraint of quark flavor 
conservation. In the microcanonical transport approach discussed here, the relations 
(^) come about due to the strong annihilation of antibaryons with baryons and the 
backward flavor rearrangement channels by detailed balance, i.e. by purely hadronic 
reaction channels in the expansion phase of the system. 



It should be pointed out, that the canonical statistical approach of Ref. 
- involving an additional parameter of dimension fm^ - well describes the strange 
and multi-strange baryon and antibaryon abundancies in Pb + Pb collisions at SPS 
energies. Thus the concept of chemical equilibrium also on the strangeness sector 
S= ±1, ±2, ±3 appears compatible with the experimental observations. However, 
it is argued here that the relations (^) are a consequence of the strong flavor 
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rearrangement reactions between formed hadrons and do not signal the presence of 
a QGP state. 

Assuming chemical equilibrium also for baryons and antibaryons with strangeness 
then (by Eq. (^9])) the A multiplicity is related to the antiproton multiplicity as 



^ ~ X - P- (51) 




All quantities on the r.h.s. of ( |FI| ) are known from the transport calculation, how- 
ever, only the interacting number of protons have to be counted in this case since 
proton spectators have to be excluded in this balance. A rather save way is to 
take into account only particles at midrapidity for |Ay| < 1, which then excludes 
spectator baryons. Thus counting only particles at midrapidity for the ratio in the 
brackets in (|51|) the A abundancy is entirely determined by particle ratios that have 
sufficient statistics in the transport calculation. The resulting {p + Y)/ < tt > ra- 



tio for Pb + Pb at 160 A-GeV is shown in Fig. |10] by the open circles with error 
bars that are due to particle statistics in the transport calculation. In this case 
the experimental multiplicity is reproduced rather well from peripheral to central 
collisions (within the error bars) suggesting the ratio A/p ^ 1 for a wide range of 
impact parameters. However, a precise experimental separation of antiprotons from 
antihyperons will be necessary to clarify the present ambiguities. 



4.2 AGS energies 

The experimental information on antibaryon production at AGS energies is rather 
scarce and limited to specific rapidity intervals. Since antibaryon yields are also 
reduced substantially as compared to SPS energies this imposes severe constraints 
on the statistics in nonperturbative transport calculations. Thus, before addressing 
any comparison to experimental data, we show in Fig. |TT] the density of 'formed' 
baryons in a central expanding cylinder of transverse radius Rt = 5 fm for Au + Au 
at 11 A-GeV in comparison to the net quark density pq/3 (dashed line). When 
comparing to Fig. ^ for Pb + Pb at 160 A-GeV roughly the same maximum net 
quark density is found for \Ay\ < 1, however, the nonhadronic phase characterized 

by 13), 

I.e. Tnonhad. 

K, 6.2 fm/c, lasts much longer due to the lower Lorentz 7- 
factor ~ 2.6. Furthermore, the density of 'formed' baryons is lower than at 
SPS energies since most of these hadronic states rescatter again on baryons in the 
medium since the formation time Tp = 0.8 fm/c is small compared to the reaction 
time roughly given by 2R/'ycm, where R denotes the radius of the Au nucleus. Thus 
'formed' baryons are reexcited to strings for a couple times during the nucleus- 
nucleus collision. This phenomenon has been addressed as 'string matter' in Ref. 



| 75[| and should not be interpreted as a state of quark-gluon plasma (QGP). 

A comparison of the HSD transport calculations with the experimental data of 
the E866 Collaboration at midrapidity p6[ is presented in Fig. O as a function of the 



participating protons Npp for Au + Au at 11.6 A-GeV/c. Whereas the proton to 
ratio is rather well described as a function of centrality (lower left part) the / 
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and K~ /'K'^ ratios are underestimated systematically for all centralities when dis- 
carding self energies for the strange hadrons (cf. Refs. |6^, |76[). The calculated 

/ ratio from the HSD calculation is 5±0.3 for all centralities rather well in line 
with the experimental observation at midrapidity. Since strangeness conservation is 
exactly fulfilled in the calculations this demonstrates that the net production of ss 
quarks by 'hard' baryon-baryon, baryon-meson and meson-meson reactions is un- 
derestimated in the transport approach as discussed in more detail in Refs. |65| . 
On the other hand the antiproton to Tr"*" ratio (lower right part) is compatible with 
the experimental ratios within the error bars indicating an approximately constant 
value of ~ 2.5 — 3 x 10~^. This roughly constant p/vr"*" ratio is a consequence of an 
approximate chemical equilibration as demonstrated in Fig. |1^ for a central collision 
of Au + Au at 11.6 A-GeV/c. Here the solid histogram corresponds to the annihi- 
lation rate of antibaryons whereas the dashed histogram stands for the backward 3 
meson fuse rate. Though the statistics are very limited in this case a net absorption 
of antibaryons, i.e. the difference in the time integrals of the annihilation rate and 
3 meson production rate, is still present in the calculations which amounts to about 
20% for central reactions and about 30% for very peripheral reactions of the total 
number of antibaryons produced in baryon-baryon or meson-baryon reactions. This 
relative net absorption, however, can only be extracted from transport calculations 
and is not a measurable quantity experimentally. 

The E877 Collaboration, furthermore, has observed a sizeable anti-flow of p's 
in Au + Au reactions at 11.6 A-GeV/c ||2^ which either indicates a strong absorp- 
tion of antiprotons on baryons or the current of comoving mesons, that fuse to BB 
pairs, and are anticorrelated to the proton current themselves. The present trans- 
port calculations reproduce these correlations, however, suffer from large statistical 
error bars (similar in size to those of the data) such that an explicit comparison is 
discarded in this work. 

As mentioned in the introduction, a high ratio of A to p of 3.6li;8 has been 
reported by the E917 Collaboration p8| for central collisions of Au + Au at 11.7 
A-GeV/c that is not understood so far. These data, furthermore, are in a qualita- 
tive agreement with the measurements from the E864/E878 Collaboration |7^. To 
estimate the A/p ratio within the present transport approach as a function of the 
centrality of the collision we employ again the relation (|51|) and invoke the parti- 
cle multiplicities at midrapidity |A?/cm| < 1- The results of these calculations are 
displayed in Fig. |T^ for Au + Au at 11.6 A-GeV/c as a function of the number 
of participating protons indicating a steady rise of the ratio with centrality. The 
hatched area in Fig. 14 demonstrates the uncertainty in the A/p ratio due to the 
limited statistics. However, the calculations for the most central collisions do not 
suggest ratios above 1.4 which is still slightly out of the range of the number quoted 
in Ref. p8| of 3.6l^;8. Thus either new theoretical concepts and/or much refined 



data are necessary to unravel this puzzle. 
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5 Summary 



In this work the conventional transport approach for two-body induced reactions has 
been extended on the formal level to n-body n ^ m reaction channels employing the 
principles of detailed balance. As a specific example of current interest the baryon- 
antibaryon annihilation problem in relativistic nucleus-nucleus collisions has been 
addressed at AGS and SPS energies where the meson densities are comparable (at 
the AGS) or much larger than the baryon densities (at the SPS) such that multiple 



meson fusion reactions as suggested in Refs. |2^, ^ become very probable. 

In order to employ the many-body detailed balance relations (as addressed in 
Section 2) a simple phase-space model for antibaryon annihilation has been presented 
that is based on flavor rearrangement channels to pseudoscalar and vector mesons 
(cf. Fig. |I]). Furthermore, a suitable covariant scheme for the calculation of such 
multi-particle reactions has been presented in Section 3 which is a straight forward 



extension of the concept proposed by Lang et al. in Ref. The method and 

its implementation in the HSD transport approach Q has been tested for a 
homogenuous system of nonstrange hadrons in thermal and chemical equilibrium 
(cf. Fig. I). 

Actual transport calculations have been performed for Ph + Ph at 160 A-GeV 
and Au + Au at 11.6 A-GeV/c, i.e. the most prominant reactions at the SPS and the 
AGS, where a couple of experimental data are available to control the dynamics. 
It is found that at both energies the meson fusion reactions are by far the most 
dominant production channel for the final antibaryons (seen experimentally) and 
that the antibaryons and baryons - at least at midrapidity - come close to local 
chemical equilibrium with the mesons. As a consequence the p/ir ratio is practically 
independent on the centrality of the collision, i.e. as a function of Apart, in line with 
the experimental observation at both energies. 

On the other hand, the approximate chemical equilibration allows to perform 
rather reliable extrapolations for the strange and multistrange antibaryon abundan- 
cies on the basis of particle ratios (cf. (^)), that follow from chemical equilibrium 
for the mesons and antibaryons. The latter ratios can be calculated with better 
statistics in the nonperturbative approach than the direct strange antibaryon abun- 
dancies. Here a roughly constant A/p ratio of ~ 1 is predicted for semi-peripheral to 
central collisions of Ph + Ph at the SPS that will be controlled soon by experimental 
data from the NA49 Collaboration. Moreover, a separation of antiprotons from an- 
tihyperons will also allow to investigate the question of antiproton self energies that 
enhance the p/vr ratio with centrality when adopting attractive scalar potentials in 



line with the analysis performed at SIS energies of ~ 2 A-GeV [T^. At AGS energies 
of 11.6 A-GeV/c the situation is less clear due to the rather low statistics for an- 
tibaryons, both theoretically and experimentally. Recall that the p/vr"*" ratio is only 
about 2.5 — 3 X 10~^ at this energy. The HSD transport calculations for the most 
central collisions of Au + Au at the AGS - assuming ( ^9|) to hold - give an upper 
limit of 1.4 for the A/p ratio, which is still slightly out of the range of the number 
quoted by the E917 Collaboration p8| of 3.61^ ;§ or related results from Ref. (j77 



Thus either new theoretical concepts and/or much refined data will be necessary to 
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unravel the antihyperon puzzle at AGS energies. 

The author likes to acknowledge continuous and valuable discussions with C. 
Greiner. Furthermore, he is indepted to E. L. Bratkovskaya and C. Greiner for 
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Figure 1: Illustration of the flavor rearrangement model for BB annihilation to 
3 mesons and vice versa. The mesons may be either pseudo-scalar or vector 
mesons, respectively. 
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Figure 2: The distribution in the final number of pions P{N.,^) for pp annihilation 
at invariant energies 2.3 GeV < ^/s < 4 GeV (short lines). The solid line is the 
gaussian parametrization ( pSj) that is fitted to the experimental data. 
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Figure 3: The annihilation cross section o"a„„. for the pp reaction as a function of 
the laboratory momentum piab from Ref. |]63| in comparison to the approximation 
(H) (solid hue). 
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Figure 4: The number of A^A^ — ppvr reactions as a function of the invariant energy 
^/s for a system in thermal and chemical equilibrium at temperature T = 170 MeV 
and fig = 0. The solid line denotes the differential number in the backward (ppvr) 
collisions, respectively. 
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Figure 5: The calculated average number of charged pions < tt > divided by Apart 
as a function of centrality for Pb + Pb at 160 A-GeV in comparison to the data from 
Ref. 



30 



0.014 
0.012 
0.010 

A 

V 0.008 

A 

V 0.006 
0.004 
0.002 
0.000 



1 ' 1 ' — 

Pb+Pb, 160AGeV 



\ - ^ no annihilation 







with annihilation + meson fusion 



with annihilation 



100 



I 

200 

part 



300 



400 



Figure 6: The ratio of antiprotons to the average number of charged pions < vr > as 
a function of centrahty for Pb + Pb at 160 A-GeV within different approximations. 
The dashed hne reflects calculations without annihilation of antibaryons, the dotted 
line includes the annihilation channels while the solid line stands for the calculations 
with both, the annihilation and meson fusion channels. 
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Figure 7: The annihilation rate BB — > Smesons (dashed histogram) for a central 
Pb + Pb collision at 160 A-GeV as a function of time in comparison to the backward 
reaction rate (solid histogram) within the HSD transport approach. 
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Figure 8: The calculated j-n (dashed), K~ jn (dotted) and p/vr (solid) ratio for 
Ph + Ph at 160 A-GeV as a function of the number of participating nucleons Apart- 
The experimental data are taken from Ref. pO[; the latter data for antiprotons 



include also the feeddown from A and respectively. The stars for Apart = 2 
denote the transport results for / vr, / it andp/ir ratios in case of A^A^ reactions 
at Tiab = 160 GeV. 
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Figure 9: The baryon density in a central expanding cylinder (see text) for | A?/|cm < 
1 in case of a central collision of Pb + Pb at 160 A-GeV in the HSD approach as a 
function of time. The solid line shows the density of 'formed' baryonic states while 
the dashed line stands for the net quark density Pg/3, while merges with the baryon 
density in the expansion phase. The dotted lines at po and 2po are drawn to guide 
the eye. 
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Figure 10: The antiproton to charged pion ratio (see text) as a function of the 
number of participating nucleons Apart for + at 160 A-GeV. The dotted hne 
denoted by f/ = corresponds to the cascade result (cf. Fig. H) while the solid line 
denoted by f/ = —100 MeV corresponds to a calculation with the attractive scalar 
antiproton potential (|48D of -100 MeV at baryon density po- The full triangles 
represent the data from Ref. [^] while the open circles correspond to a cascade 
calculation including the feeddown from antihyperons according to Eq. (^TJ). The 
errorbars are due to the limited statistics in the transport calculation. The star for 
A-part = 2 denotes p/vr ratio for A^A^ collisions at Tiah = 160 GeV. 
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Figure 11: The baryon density in a central expanding cylinder (see text) for | A?/|cm < 
1 in case of a central collision of Au + Au at 11.6 A-GeV/c in the HSD approach 
as a function of time. The solid line shows the density of 'formed' baryonic states 
while the dashed line stands for the net quark density Pg/3, while merges with the 
baryon density in the expansion phase. 
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Figure 12: The /tt~^, K~ /tt~^, p/7c~^ and p/n'^ ratio at midrapidity for Au + Au 
at 11.6 A-GeV/c as a function of the number of participating protons Npp. The sohd 
hne correspond to the HSD transport calculation in the cascade mode for mesons 
and antibaryons (open circles) while the experimental data (full symbols) are taken 
from p6| . 
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Figure 13: The annihilation rate BB Smesons (dashed histogram) for a central 
Au + Au collision at 11.6 A-GeV/c as a function of time in comparison to the 
backward reaction rate (solid histogram) within the HSD transport approach. 
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Figure 14: The K/p ratio as a function of centrality for Au + Au at 11.6 A-GeV/c 
in the HSD approach within the approximation (^l]). The shaded area corresponds 
to the uncertainty in the statistics of the transport calculations for the different 
particle abundancies at midrapidity. 
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